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Abstract 

Critical phenomena in the two-dimensional Ising model with a defect line are 
studied using boundary conformal field theory on the c = 1 orbifold. Novel 
features of the boundary states arising from the orbifold structure, including 
continuously varying boundary critical exponents, are elucidated. New fea- 
tures of the Ising defect problem are obtained including a novel universality 
class of defect lines and the universal boundary to bulk crossover of the spin 
correlation function. 
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Boundary conformal field theory |l| is of considerable current interest in string theory, 
classical two-dimensional critical phenomena and quantum impurity problems. || In the case 
of the c = 1 free boson, compactified onto a circle or radius r, Dirchlet (D) and Neumann 
(N) boundary conditions, in which the field or its dual obey tp = y? or <P — <Poi have been 
well studied. 0] By symmetry, the dimensions of boundary operators are independent of 
<Po, <pQ. Here we extend this analysis to the c = 1 Z 2 orbifold || in which (p and —ip are 
identified, showing that now the boundary dimensions depend on these parameters and that 
special features occur at the orbifold fixed points, (p = 0, nr, <p = 0, Tr/2r. These results 
are used to understand the critical behavior of a defect line in the Ising model. 

We consider a two-dimensional Ising model at its critical temperature, with a (horizontal) 
defect line. Modifying the vertical couplings across the defect line, for a square lattice on a 
torus, yields the classical Hamiltonian: 

£ = ~~ X] [Jl <T i,j <T i+l,j + h&ijffij+l] 

i=l,M j=l,N-l 

- X! [Jl a i,N^i+l,N + J&i,NCTi,l]- (1) 

i=l,M 

Various exact results have been obtained on this model. P,[7| When the bulk couplings, 
Ji and J 2 are tuned to the critical point, the defect model exhibits a line of critical points 
depending on the details of the defect, i.e. the value of J in Eq. ([[]). The spin-spin correlation 
function decays along the defect line with a critical exponent which varies continuously along 
the line of critical points. Other types of defect lines are possible, including modification of 
horizontal couplings, introducing multi-spin couplings, imposing magnetic fields at the line, 
etc. 8 other universality classes can be readily obtained by cutting the system at the defect 
line and imposing independent boundary conditions, spin- up, spin-down or free on the 2 
sides. 

To apply boundary conformal field theory (CFT) to this problem, we must "fold" the 
system at the defect line, || (See Fig. 1) obtaining 2 Ising variables at each point, coupled 
only at the boundary, a special case of the Ashkin- Teller (AT) model. The corresponding 
CFT is hence the c = 1 Z 2 orbifold. || We show that all known universality classes of 
boundary conditions in the Ising defect line problem, and some additional ones not previously 
known, correspond to D and N boundary conditions on the orbifold model. The universal 
groundstate degeneracy, g, is calculated for all boundary conditions and the relative 
renormalization group stability is discussed. We have calculated the spin correlation function 
with the spins at arbitrary points relative to the defect line and each other for arbitrary 
boundary conditions; we give the result for D boundary conditions here. Details and further 
results will be given in a longer paper. [[UJ 

The Hilbert Space for the AT model with periodic boundary conditions contains two 
sectors, in which the free boson obeys either periodic or twisted boundary conditions. || 
Boundary states are constructed, in general, using both sectors. (Related discussions in 
string theory context are found in WM.) For a periodic boson, boundary states are con- 
structed using the oscillator creation operators for right and left movers of momentum 
±27in/j3, a) Rn and a} Ln where (3 is the length of the periodic spatial interval. In addition the 
zero-mode states, \(w, k)) are used. Here w and k are integers which label winding numbers 
in the space and time direction respectively. The left and right-moving energies of these 
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states are (2tt/(3)x and (2ir//3)x with x = (rw + k/2r) 2 , x = (rw — k/2r) 2 . The Dirichlet 
boundary state which is an eigenstate of f(u) with eigenvalue ipo is: [|J 



|/%o)> = ^exp[-£aL4j E e-*">/'|(0,*)>- (2) 

V ±T n =l fc=-oo 

Orbifold model boundary states must be invariant under ip — > —ip; hence it is necessary to 
combine the D states with eigenvalue <p and — p : 

\D o (<p )) = ±[\D(<p )) + \D{-<p ))], (3) 

It is useful to consider a torus of circumferences (3 and 21 with two D defect lines at dia- 
metrically opposite locations, separated by I. After folding, this corresponds to the finite 
cylinder geometry considered in Ref. 0. The partition function, at inverse temperature (3 
can be written in terms of the boundary states as: 



W = (DoMle^fi \DoM) = Z(p - <f/ ) + Z(p + tp' ), (4) 



where: 



Z{y Q -y' Q ) = {D{y Q )\e- lH e\D{y' Q )) 



1 oo 

with 77(g) = q ,1 ^ 24 n^i(l — (f 1 ) and q = e~ n ^' 1 . Z(ip — ip' ) is the partition function for a 
periodic boson. Note that all states in the Boltzmann sum have integer multiplicity and 
that, in the case ip = <p' Q , the groundstate occurs with multiplicity one. It is this physical 
condition which fixed the normalization of the boundary state in Eq. fl2|) and In the 
case of equal boundary boundary conditions, the finite-size energies give the dimensions of 
boundary operators fj]. For a periodic boson the boundary operator spectrum can be read 
off from Eq. (|^): x — 2r 2 n 2 + integer. It is independent of (po, as it must be by symmetry. 
However, this is not the case for the orbifold boson; instead, from Eq. @, the boundary 
scaling dimensions vary continuously with ip , corresponding to a line of fixed points. 

Also note that for the orbifold case when <p is at one of the fixed points, or nr, the 
groundstate occurs with multiplicity two, signalling that this is not a well-defined boundary 
state. To obtain well-defined Dirichlet boundary states in these two cases we must add a 
component from the twisted sector. Imposing the twisted boundary condition ip(a + 0) — 
—if (a), leads to oscillators with half- integer momenta, a n+ \/2- There are no zero-modes 
in this case and the constant part of f must have the value or 7rr, corresponding to two 
different oscillator groundstates, |0)r and |7rr)r. Thus we can only construct Dirichlet states 
in the twisted sector for these values of y?o- These states are: 

00 

\D(<P0)t) = eX P [- £ 4(n+l/2) a fl(n+l/2)] \<Po)t, (6) 
n=0 

where <y9 = or 7rr only. There are four possible D boundary states in these cases, which 
give integer multiplicities when combined with the other D-states and give unit multiplicity 
for the groundstate when combined with themselves. These are: 
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\DoM±) = 2- 1/2 | J D(^o)) ± 2- 1 / i \D( Vo ) T } } (7) 

for ifo = 0, ttt only. 

In a similar manner we can construct N-states. For a periodic boson these are: 



\N{^)) = V~r £ e-^expt+E^UUl^^)). (8) 

71=1 



Note that in this case we use the other type of zero-modes which vary in the space-direction, 
and that the sign in the exponential for the oscillator factor is reversed. |iV(<^o)) is an 
eigenstate of the dual field, <p(cr), of radius l/2r, with eigenvalue <pq. We may regard <p as 
an orbifold variable with fixed points at and tt/2t. Thus the orbifold N states are: 

\N o (0 o )) = ^[\N(<p o ) + \N(-<p o ))}, (9) 

for < (fio < it /2r. We can construct 4 special N-states for <^ = 0, 7r/2r by analogy 
with the D construction, using the twisted sector. These are constructed from the linear 
combinations of oscillator groundstates (\0)t ± l 7rr )T), respectively. 

These 2 continuous lines, < (p < ttt and < y5 < 7r/2r plus 8 the discrete points 
of boundary states obey the Ishibashi condition [T(er) — T(er)]|A) = where T and T are 
the left and right-moving components of the energy-momentum tensor. All partition func- 
tions constructed from any pair of these boundary states contain only non-negative integer 
multiplicities. The states constructed from the same boundary state at both boundaries all 
have unit multipicity for the identity operator. We conjecture that for generic values of r 
this is the most general set of boundary states satisfying these conditions. The analogous 
statement has been proven for a periodic boson. [12[ 



We now apply the r = 1 case of this general analysis of orbifold boundary CFT, in to 
the Ising model with a defect line. The D-boundary conditions correspond to the integrable 
defect line of Eq. ([I]). By comparing the exact spectrum J7| with that of the CFT we obtain: 

tan(^ - tt/4) = Sinh[(1 ~ J/J2)J2/r] (0 < <p < vr), (10) 
^° ' ' sinh[(l + J I J 2 )J 2 /T] y ^ h y > 

while sinh ( J\/T) sinh(J 2 /T) = 1 so that the bulk Ising model is critical. As we vary J ' j J 2 
from — oo to oo, ip decreases monotonically from ir/4 — tan _1 (exp(2J2/T)) (between 2>tt/A 
and 7r) to 7r/4 — tan _1 (exp(— 2J 2 /T)) (between and tt/4). Only in the extreme anisotropic 
limit J 2 /J\ — > 0, does ip approach the endpoint values tt and 0. 

There are additional discrete universality classes of defect lines in which different bound- 
ary conditions are imposed independently on the two sides of the line. Cardy has shown [l 
that there are only three universality classes of boundary conditions for an Ising model on a 
semi-infinite plane: spin-up (j), spin-down (j) and free (/). We denote the boundary state 
corresponding to a f boundary condition on the left of the defect and a | on the right as 
| t|), etc. The \ff) state simply corresponds to the Dirchlet state with J / J 2 = 0, y? = tt/2. 
By direct comparison of the partition functions, we have verified that the four boundary 
states corresponding to \ and { boundary conditions are: 
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I TT) = \D o (0)+), | IT) = \D (ir)+), (11) 

with the spin-reversed states given by the corresponding "-" orbifold boundary states. It is 
interesting to note that the unphysical boundary state obtained by extrapolating ip — > 
without including the twisted sector can be written: 

lim \D M) = \D o (0)+) + |A>(0)-) = | TT) + | II). (12) 

This limit corresponds to an infinitely strong coupling across the defect line and also a limit 
where the other vertical couplings J%jT — > while the horizontal couplings, J\/T — > oo. 
Thus it seems reasonable to suppose that the two horizontal chains of spins which are coupled 
across the defect line, (i, N) and (i, 1) get locked into a perfectly ferromagnetically aligned 
state corresponding to | TT) + I ID- Similarly, the limit (p — > n gives the perfectly aligned 
antiferromagnetic state. 

The other four possibilities for imposing different boundary conditions on the two sides 
of the defect line, |/ f), etc correspond to the four endpoint N boundary conditions in a 
similar way. The nature of the N-line of boundary states is commented on below. 

We now turn to the RG stability of these various universality classes of boundary con- 
ditions. This can be addressed by checking for relevant boundary operators allowed by 
symmetry, and also by invoking the ^-theorem which states that the groundstate degen- 
eracy, g, always decreases in flows from less stable to more stable fixed points. We first 
consider the continuous line of D-fixed points, with < ip < ir. This has the value, g = 1, 
along the entire line, indicating that this is indeed a line of stable fixed points; no flow can 
occur along the line. We can read off the dimensions of all boundary operators from the 
partition function of Eq. The only non-trivial relevant boundary operators (of dimen- 
sion x < 1) have dimensions 2((po/ir) 2 and 2((Pq/tt — l) 2 . They can be shown to correspond 
to <7i ± 0^, the sum and difference of the Ising spin operators on the two sides of the defect 
line. If we consider a defect line with the Ising Z 2 symmetry, then the relevant boundary 
operators are prevented from occuring by symmetry. Thus the D-boundary conditions are a 
stable line of fixed points which should attract generic defect lines with Z 2 symmetry. That 
is, we may consider arbitrary vertical and horizontal couplings near the defect line (which 
preserve horizontal translational symmetry) and always expect to renormalize to a D fixed 
point with some value of ipo. On the other hand, if we break the Z 2 symmetry, for example 
by appling a magnetic field along the defect line, then we destablize these fixed points. It is 
natural to expect that we then renormalize to one of the endpoint D-boundary conditions 
corresponding to independent spin up or down boundary conditions on the two sides of the 
defect line. Taking over the results for a semi-infinite Ising system, |1],|9| we conclude that 
these boundary conditions have g = 1/2, and no relevant operators. They are the most 
stable fixed points. 

The general N-states have g = y/2 and have relevant operators, of dimension x = 1/2 
which correspond to a product of the Ising spin operators on the two sides of the defect 
line. Thus, they represent unstable critical points which could only be achieved by imposing 
additional symmetries at the defect line. |L0| We expect that adding other couplings would 



produce a flow away from the N boundary conditions to D boundary conditions (g = 1), if 
the Z 2 symmetry is preserved, or to fixed spin boundary conditions (g = 1/2) if not. 



5 



Finally, we turn to the critical spin-spin correlation function for an infinite system with 
a defect line of D-type, for arbitrary strength of the defect. We introduce a complex spatial 
co-ordinate for the folded system, z = x + iy, where x, y are the two spatial co-ordinates; the 
defect is at y = so y > after folding. (If the Ising model does not have Ji = J 2 , square 
symmetry, it is neccessary to rescale one of the co-ordinates.) The folding introduces two 
different spin operators, <Ti and <ji at each point. This universal bulk to boundary crossover 
function, depending in a non-trivial universal way on the cross-ratio: x = | {z\— Z2) / {z\— Z2) | 2 , 
can be determined by standard methods in terms of the boundary state. 0,13]. Using the 



explicit construction of the AT model conformal blocks by Zamolodchikov [J14|| , and the 
D-boundary states given above, we obtain: 

/ l \ V8 



<ff ' ( ^ fe)> = te] M^m Me2 ^ u{x)) ' (13) 

for two spins on the same side of the defect line (before folding). Here dj are elliptic theta 
functions (as defined in Ref. ||) and u(x) is defined by x = [^(u) / '$ 3 (u)] 4 . When the two 
points are on opposite sides of the defect line before folding we obtain: 



Using Mathematica, we have plotted, in Fig. (2), the correlation functions for the two 
points at unit distance from the boundary and separated by an arbitrary distance from each 
other. In the short-distance limit, (aia\) converges to a unique power-law, 1/4, which is 
independent of the defect strength. This is the expected bulk limit. In the large-distance 
limit, the correlation function is governed by another exponent, (2<^o/tt) 2 , which depends 
on the defect strength. This is the boundary limit. Our result interpolates between these 
two limits. In the "short-distance" limit, (<7i0"2) converges to a constant which depends on 
the defect strength. This actually corresponds to two spins located symmetrically about the 
defect line. In general, (ctio^) is smaller than {<3\<Ji) for the same defect strength and the 
same (horizontal) distance, as expected. Nevertheless, they asymptotically converge to the 
same power-law function in the large-distance limit with the constant prefactor for J 7^ 0. 

We thank F. Lesage for useful information. This work is partly supported by NSERC. 
M. O. thanks the Killam memorial fellowship for financial support. After most of this work 
was completed, we were informed that A.W.W. Ludwig has obtained some related results 
independently. 
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FIGURES 





Defect Line 

FIG. 1. The folding of the Ising model on a cylinder to a c = 1 theory on a strip. We fold at 
the defect line and also at the line on the opposite side. These lines correspond to the boundary 
in the folded system. 
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(b) Correlation function of two spins in different layers 
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FIG. 2. The two-spin correlation for various strength of the defect. We show the result for 
(a) (aiai) and (b) {^1^2) for ipo = (strong coupling and anisotropic limit), 0.17T, 0.2-7T, 0.25-7T 
(no defect) , 0.37T, 0.47T and 0.57T (free boundary condition). They are shown as a function of the 
(horizontal) distance r, in a log-log plot. 
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